To APPEAR IN "XIII lAU LARIM (2011) 



RevMexAA(SC) 



THE DARK AND STELLAR MASS ASSEMBLY OF GALAXIES 
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ABSTRACT 

The emerging empirical picture of galaxy stellar mass (Mg) assembly shows that galaxy population buildup 
proceeds from top to down in M^- By connecting galaxies to ACDM halos and their histories, individual 
(average) Ms growth tracks can be inferred. These tracks show that massive galaxies assembled their Mg the 
earlier the more massive the halo, and that less massive galaxies arc yet actively growing in A/g, the more active 
the less massive is the halo. The predicted star formation rates as a function of mass and the downsizing of the 
typical mass that separate active galaxies from the passive ones agree with direct observational determinations. 
This implies that the ACDM scenario is consistent with these observations. The challenge is now to understand 
the baryonic physics that drives the significant and systematical shift of the stellar mass assembly of galaxies 
from the mass assembly of their corresponding halos (from halo upsizing to galaxy downsizing). 
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1. INTRODUCTION 

Galaxies are fascinating astronomical systems. 
On the one hand, the complexity of stellar popu- 
lations, ISM, AGN, and dark matter is integrated 
inside them. On the other hand, galaxies are the 
lost link between the early universe and the observed 
world of astronomical objects; as such they trace the 
way cosmic structure emerged in the universe and are 
the basis of unique methods for determining several 
cosmological parameters. One of the current chal- 
lenges in the study of galaxies is determining and un- 
derstanding how did they assemble their stellar and 
dark halo masses. The main observational imprint 
of the complex galaxy assembling processes is prob- 
ably the stellar-to-dark mass ratio, Fs = Mg/Mh, as 
a function of mass and redshift z. 

The growth of Mg can be due (1) to local 
star formation (SF) fed by available or infalling 
gas and/or (2) to capture of other stellar sys- 
tems. In the former case, SF is driven by in- 
ternal physics, mainly b y self-regulation processes 
in the galaxy disk (e.g., Firmani fc Tutukovl Il992 : 



Firmani fc Avila-Reesel2000l) . or is induced by exter 
nal interactions and mergers -mainly those with high 
gas content. In the latter case, Afg can grow by accre- 
tion of small satellites and tidal debris or by major 
mergers. Even more, the mentioned processes imply 
feedback effects, like SN-driven outflows and AGN- 
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driven intrahalo gas reheating/ejection effects that 
inhibit further Mg growth. Which of these channels 
and their respective feedback effects did dominate in 
the Ms growth of galaxies as a function of their halo 
masses, types, and environment? How does proceed 
the stellar and dark mass assembly of galaxies ? 

2. DIFFERENT APPROACHES TO GALAXY 
FORMATION AND EVOLUTION 
2.1. The inductive or 'archaeological' method 
The approaches for studying galaxy formation 
and evolution have changed over time. In the past, 
the main approach was the 'archaeological' method: 
ages and star formation rate histories of our and 
other local galaxies arc reconstructed by means of 
stellar population synthesis and chemical evolution 
models from the observed spectro-photometric and 
chemical properties. The application of this method 
to modern data and large galaxy surveys have al- 
lowed to obtain key constraints to galaxy mass as- 
sembly as a function of mass, galaxy type, and en- 
vironment (e.g., Thomas et al. 2005; Panter et al. 
2007; Asari et al. 2007). 

2.2. The deductive or ab initio method 
With the consolidation of the popular A Cold 
Dark Matter (ACDM) cosmology in the last twenty 
years, a powerful theoretical background for galaxy 
formation and evolution appeared. According to 
this framework, cosmic structures emerged from the 
gravitational evolution of CDM-dominated primor- 
dial perturbations. The distributions, inner prop- 
erties, and evolution of the collapsed CDM struc- 
tures (halos) were predicted in detail by cosmological 
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N-body simulations. Ab initio (deductive) models, 
as well as full numerical simulations including hy- 
drodynamics, were developed in order to follow the 
complexity of baryon physics during the process of 
galaxy formation and evolution inside the growing 
CDM halos. Several of the physical ingredients and 
parameters of the models, as well a some of the sub- 
grid parameters of the simulations, are actually de- 
termined/calibrated from mainly local observations. 
The deductive approach has enormously contributed 
to our understanding of galaxy formation and evolu- 
tion in the cosmological context though many ques- 
tions remain yet unsolved, including the validity of 
the backbone of this approach itself: the existence of 
the elusive CDM. 

The halo mass growth in the ACDM cosmology 
is on average hierarchical, from bottom to up. How 
is the Ms assembly of galaxies formed inside the 
CDM halos and how does it compare with obser- 
vations? A sharp test for the ACDM scenario may 
arise from this comparison. For this, however, sev- 
eral theoretical and observational "biases" should be 
first well understood. Neither the galaxies trace di- 
rectly the assembly of their halos -here is implied 
all the complex baryon physics-, nor the observed 
galaxy populations at different redshifts offer a di- 
rect way to determine the individual evolution of 
galaxies -instead they are inferred and the results 
are strongly dependent on selection effects and bi- 
ases. It should be said that misunderstandings and 
oversimplifications of the scenario have led to some 
incorrect interpretations. For example, it is com- 
mon to hear that since in the ACDM hierarchical 
scenario dark halos assemble through violent ma- 
jor mergers ("walnut tree"), therefore galaxy forma- 
tion is merger-dominated, something that could be 
in conflict with the observed large abundance of spi- 
ral galaxies or the non-negligible fraction of bulge- 
less galaxies. In fact, most of the mass of galaxy- 
sized ACDM halos (and most of halos) was (were) 
not assembled by major mergers but by smooth ac- 
cretion and minor mergers ("pine tree"; Maulbetsch 
ct al. 2007; Genel et al. 2010). Besides, only a 
fraction of halo-halo mergers end producing galaxy- 
galaxy mergers. By analyzing the results of the SAM 
implemented in the Millenium simulation (De Lucia 
et al. 2005), in De Rossi et al. (2009) it was con- 
cluded that only « 12% of Milky Way- type galaxies 
experienced a major merger during their lifetimes. 

2.3. The empirical approach 

Thanks to the vertiginous advance of instrumen- 
tal facilities and multi-wavelength large-area survey 



programs, the empirical or "look-back time" ap- 
proach has come to stay. The main photometric and 
spectral properties of whole galaxy populations ob- 
tained at different z allow us to infer Mg, SF rate 
(SFR), and other physical properties of galaxiefH. 
Given the large numbers and high sensitivities of 
current surveys, relatively complete samples char- 
acterized by the mentioned above properties down 
to nearly dwarf galaxy masses and out to high red- 
shifts become possible (e.g., Ms>10® at z « or 
Ms > 3 10^ at z ~ 2). Thus, progress has been made 
in determining locally and at higher z: (1) the total 
galaxy stellar mass function, GSMF, and its dissec- 
tion into blue/star- forming and red/passive compo- 
nents, (2) the SFRs of galaxies as a function of their 
mass traced by different indicators, from the UV to 
the radio, and (3) the galaxy rate of merging at dif- 
ferent masses and epochs, although with still very 
large systematic uncertainties. In some works, these 
determinations have been obtained as a function of 
environment. 

Based on all these observational data, a purely 
empirical picture of Ms assembly as a function of 
mass, type and environment is emerging now. A 
general result is that ^ 50% of the local stellar mass 
density was assembled since z = 1, and ^ 90% since 
z = 3.5 (e.g., Perez-Gonzalez et al. 2008). Possibly 
the main new concept distilled from the empirical 
pic ture is that o f cosm ic downsizing, a term coined 
by ICowie et al.l (|l996l ) to describe the decline with 
cosmic time of the maximum rest-frame iiT— band 
luminosity of galaxies undergoing active SF. This 
term has been more recently used to describe a num- 
ber of trends of the galaxy population as a function 
of mass that imply in general galaxy buildup from 
top to down. However, these different trends are 
actually related to different astrop hysical phenom- 



ena a nd galaxy evolutionary stages (jFontanot et al, 
l2009l) . From the most general point of view, the 
many downsizing manifestations can be separated 
into those that refer to the evolution of: 

• (A) Massive galaxies, which today are on av- 
erage red and passive (quenched SF). Observa- 
tions show that the high-mass end of the GSM F 
was in place since high z (e.g., Fontana et al. 
2004,2006; Drory et al. 2005; Cimatti et al. 
2006; Marchesini et al. 2009; Perez- Gonzalez 
et al. 2008), and evidence a decrease with cos- 
mic time of the characteristic mass at which the 



^For several of these inferences, stellar population synthesis 
models should be used. These models had a vigorous advance 
in the last years, however, several uncertainties remain yet 
unsolved (for a recent review see Bruzual 2010). 
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SFR is dramatically quenched or at which the 
GSMFs of early- and late-type galaxies cross, 
i.e. less and less massive migrate with time 
to the red sequence (e.g., Buudy et al. 2006; 
Borch et al 2006; Bell et al. 2007; Hopkins et 
al. 2007; Drory & Alvarez 2008; Vergani et al. 
2008; Pozzetti et al. 2010). This is in line with 
'archaeological' inferences for local galaxies that 
evidence an early and co heval mass assembly for 
massive ellipticals (e.g., Thomas et al. 2005[) . 



(B) Less massive galaxies, which are on av- 
erage blue and star-forming. Although the 
observations of less luminous galaxies suffer 
incompleteness as z increases due to the flux 
limits, at least up to z ~ 1 — 2, most of 
observational studies have found that the spe- 
cific SFR (SSFR = SFR/Ms) of galaxies with 
Ms S3 X 10^° Mq, which are mostly late- type, 
star-forming galaxies, is surprisingly high even 
at z ^ (late mass assembly) and, on average, 
the lower the mass, the highe r the SSFR 
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By combining the evolution of the GSMF with 
the measured SFR-Mg relations at different z, the 
contribution of local SF and galaxy accretion (merg- 
ers) to the stellar mass build up can be constrained 
(Bell et al. 2007; Drory & Alvarez 2008; Pozzetti 
et al. 2010). These and some direct -but yet 
limited- observational determinations of the merg- 
ing rate (e.g., Lotz et al. 2007; Bundy et al. 2009) 
show that the former channel completely dominates 
in low- and intermedium-mass galaxies at all epochs, 
while (dry) mergers may play a moderate role for 
massive (Mg^ 10^^ Mq, mainly red) galaxies at later 
epochs (z'zl)- 

Although important efforts have been made al- 
ready in order to constrain the dynamical evolution 
of galaxies (e.g., from observational studies at high 
redshifts of the Tully-Fishcr and Faber- Jackson rela- 
tions, galaxy-galaxy weak Icnsing, etc.), direct con- 
straints of the galaxy- halo connection as a function of 
mass are yet very limited. This kind of observations 
together with those regarding the gas mass of galax- 
ies at high z are necessary to complete the whole 
picture of galaxy stellar, baryonic, and dark mass 



assembly. In the meantime, semi-empirical meth- 
ods were introduced in order to get the whole evo- 
lutionary connection in the context of the ACDM 
scenario and constrain in this way the average in- 
dividual trends of the galaxy evolution process as 
a function of mass (e.g., Conroy & Wcchsler 2009; 
Hopkins et al. 2009; Firmani & Avila- Reese 2010). 

3. CONNECTING GALAXIES TO DARK HALOS 

In the semi-empirical approach, the information 
provided by direct observations at different epochs 
is combined in a statistical way with predicted prop- 
erties for the ACDM halos in order: (1) to attain 
a connection between galaxies and halos, and (2) to 
infer individual galaxy 'hybrid' evolutionary tracks 
(GHETs). An schematic idea of the approach is 
shown in Fig. [T] and explained below: 

(1) Global gal£Lxy-halo connection.- Given 
the measured cumulative GSMF and the cumulative 
halo mass function, HMF (both at the same epoch 
e.g., z = 0, black solid curves), and assuming a one- 
to-one correspondence between Mg and Mh, both 
functions are matched, $g(> Ms) = $h(> Afh), for 
finding the Mh corresponding to a given Mg (dot- 
ted arrows, for two different masses, in Fig. [IJ, i.e. 
the Mg-Mh or Fs = M^/Mh - Mh relations are in- 
ferred (solid bell-shaped curve in the inset). This 
method, called the abundance matching technique 
(AMT; e.g.. Vale & Ostriker 2004; Kravtsov et al. 
2004; Behroozi, Conroy & Wechsler 2010 -BCW10-, 
and for more references see therein), makes a mini- 
mum of assumptions and has proven to be effective 
and practical. The halo masses of galaxies can be 
actually determined by direct methods as galaxy- 
galaxy weak lensing and satellite kinematics (e.g., 
Mandelbaum et al. 2006; More et al. 2010). Never- 
theless, in current studies, due to the low signal-to- 
noise ratios, results can be obtained only by stacking 
a large number of galaxies; this introduces biases and 
significant statistical uncertainties in the inferences, 
and limits these inferences to relatively small mass 
ranges. However, in the mass ranges where compar- 
isons arc feasible, these direct methods, the modcl- 
dcpendcnt indirect methods (e.g., the Halo Occupa- 
tion Model and the Conditional Luminosity Function 
formalism), and the AMT, give local Fg-Mh rela- 
tions compatible among them within a factor of 2 
(see BCWIO; Moster et al. 2010; Guo et al. 2010; 
More et al. 2010; Rodn'guez-Puebla et al. 2011). 

The obtained stellar mass fractions, Fg, are very 
low with respect to the universal baryon fraction 
{Fh_u = Qh/riM ~ 0.16), i.e. galaxy SF inside 
ACDM halos seems to be a very inefficient process. 
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Fig. 1. Halo and galaxy stellar mass cumulative functions at z = (solid black curves). Given a Mg, its corresponding 
Afh is found by matching the mass functions (the dotted black arrows show this for two masses). A relation between 
(or Mb/ Myi) and Mh is constructed this way (solid black curve in the inset). The same exercise can be repeated at other 
redshifts (e.g., z = 1, dashed blue curves). Given the halo average MAH, the mass at z — 1, Mh(l), that had a halo of 
a given mass at z = 0, Mh(0), can be calculated (AMh). Then, by using again the AMT, the Ma(l) corresponding to 
Mh(l) is obtained (the dotted blue arrows show this for two masses), and AMs = Afs(O) — Ms{l) gives the individual 
average Afs track (GHET; the red arrows in the inset show two GHET-based tracks in the Afs/A/h-Afh diagram). 



Besides, the efficiency is strongly dependent on mass: 
it peaks around Mh = 8 10"'^^ AIq, decreasing signif- 
icantly towards lower {Fg cx My,"', a « f.25 ) and 
higher {Fg oc Mh^ & « -0.6) masses. The GSMF 
determinations at higher redshifts allow now to use 
methods hke the A MT to infer the Mg-Aih relation 
at dif ferent epochs ( Moster et al.ll2010t Wang fc Jing 
2010l BCWIO). The latter authors extended the 
AMT out to z w 4 and found that (a) the Mh at 
wich Fs peaks shifts little to higher masses (by « 0.6 
dex out to z w 4), and (b) the peak value remains 
roughly constant. For masses below (above) the peak 
at a given epoch, Fg is smaller (slightly larger) than 
at earlier epochs for a given (see Fig. O dashed 
blue curves). This implies that galaxy SF in mas- 
sive halos should have been slightly more efficient 
in the past, while in low-mass halos it should have 
been even less efficient. This is a manifestation of 
the cosmic downsizing discussed above. 

(2) Individual M^ and Ms growth tracks. - 

In the ACDM scenario, the individual mass aggre- 
gation histories (MAHs) of halos are known. By 
using them in combination with the semi-empirical 
Mg-Mh relations at different z, the corresponding 
Ms evolutive tracks (GHETs) can be inferred (red 
arrows in the inset of Fig. [1] IConrov fc Wechsler 



20091: iFirmani fc Avila-Reesell2010l hereafter FAlO). 
In FAlO, average halo MAHs and a parametrization 
of the BCWIO Ms-Mh relations up to z = 4 were 
used to infer individual average GHETs. The re- 
sults are encouraging and show that at each epoch 
there is a characteristic mass that separates galax- 
ies into two populations (Fig. [2]): (a) galaxies more 
massive than Ms{z = 0)« 3 x 10^° Mq are on aver- 
age quiescent /passive (their Ms growth slowed down 
or stopped completely), besides the more massive is 
the galaxy, the earlier it has transited from the active 
(blue, star-forming) to the passive (red, quenched) 
population ('population downsizing'); (b) galaxies 
less massive than Ms{z = 0)« 3 x 10^° Mq are on av- 
erage active (blue), and the less massive the galaxy, 
the faster its late Ms growth, driven likely by local 
SF ('downsizing in SSFR'). 

The Ms{z) average tracks (GHETs) inferred in 
FAlO are shown in Fig. |3] (solid red curves). Their 
corresponding halo MAHs arc also plotted (dot- 
dashed blue curves) but for comparative reasons, 
they were shifted vertically in such a way that each 
MAH coincides with its related GHET at z = 0. 
The shapes of the average stellar and halo mass as- 
sembling histories are quite different. For galaxies 
with Ms{z) < 3 10^° Mq, their halo MAHs at later 
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Fig. 2. The Ms/Mh-Mh functions inferred by tlie AMT at four z (dasiied blue curves; BCWIO, see FAIO for modifica- 
tions). Solid red curves are tlie average Ms growth tracks inferred in FAIO: in small halos, Ms assembles faster at late 
epochs as smaller is Mh, while in large halos, Ms stopped its growth, the earlier the more massive the halo. The two 
dotted curves are tracks corresponding to two disk galaxy evolutionary models. 

z 

epochs grow slightly slower as smaller is the mass, 
while their stellar GHETs grow much faster. For 
AIs{z~Q)> 3 10^° Mq, as the system is more mas- 
sive, the stellar assembly of the galaxy occurs earlier 
in time with respect to the corresponding halo. Only 
for systems that attained at a given epoch masses 
Ms « 3 - 5 10^° Mq (Mh « 1 - 2 1012 Mo), both 
the Ms and Mh assembly history shapes are similar; 
these systems are namely those in the peak efficiency 
(maximum Afg/Mh ratio, see Fig. [2]). 

3.1. Predictions 

SSFR histories.- By using the GHETs, the spe- 
cific Ms growth rate histories, Ms{z)/Ms{z), can be 
calculated. If working hypothesis one assumes 
that the Ms growth of a given galaxy is only due 
to local SF, then its SSFR history is Ms(^)/Ms(z) 
divided by (1 — R), where R ?» 0.4 is the gas recy- 
cling factor due to stellar mass loss. In Fig. 0] we 
reproduce the individual average SSFR tracks in the 
SSFR vs Ms diagram as obtained in FAIO (thin red 
solid lines) : the SSFR of low-mass galaxies decreases 
with time slowly, their late Mg growth being very 
fast, while for large galaxies, the SSFR decreases 
very fast and their Ms growth practically stops at 
late epochs. 

The thick solid lines correspond to isochrones 
at z = 0,1,2,3 and 4, from bottom to top, re- 
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Log(1+z) 

Fig. 3. Average Ms growth tracks (GHETs) as a func- 
tion of z (solid red curves) and their corresponding halo 
MAHs (dot-dashed purpure curves) shifted vertically in 
order Ms{z — 0)=Mi^{z = 0). Galaxies and halos assem- 
ble their masses in a very different way (from FAIO). 

spectively. These isochrones are the SSFR-Mg re- 
lations at the given z and can be compared with the 
direct determinations mentioned in §§2.3. Within 
the current large uncertainties and sample selec- 
tion effects, the GHET-based and directly deter- 
mined average SSFR-Ms relations out to z ^ 2 
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Fig. 4. Predicted average SSFR-Mg tracks (assuming 
negligible Ms growth by mergers; solid thin red curves, 
from FAIO). The thick solid curves connect the individ- 
ual tracks at a given z (isochrones). The crosses show 
when a given track strongly falls -the galaxy becomes 
passive. The mass Ms at this epoch is called Mmigr- 



are in reasonable agreement. If any, the former are 
shallower than the latter at the low-mass side; at 
low 2, the predicted SSFRs for intermediate masses 
(Ms ~ 0.7 - 2 10^° Mq) are slightly larger than the 
averages of those directly measured. It is important 
to note that the fast decreasing in the SSFR-Mg rela- 
tions (isochrones) seen at the high-mass side is asso- 
ciated to passive (red) galaxies. Most of the SSFR- 
A/g relations presented in the literature refer to only 
star-forming (blue) galaxies; for masses lower than 
A/g « 2 — 3 lO^"Af0, blue galaxies by far dominate 
in number, but at high masses passive (red) galaxies 
become dominant (e.g., Bell et al. 2003; Salim et al. 
2007; Pozzetti et al. 2010). In the works where the 
samples were partitioned by color, the SSFR-Mg re- 
lation ofthfiredgaJaxiesJsJn^ very steep 
(e.g- lBell et"al1l2007t iKarim et al.ll2010[ ). 

Downsizing of the migration (quenching) mass.- 
The GHET-based SSFRs can be used to calculate 
the z at which the SSFR function of Afg for 
a given track decreases dramatically (the galaxy mi- 
grates to the passive sequence). The crosses in Fig. |4] 
show, for each track, the Mg attained by the galaxy 
at this epoch, i.e. at each epoch there is a mass 
Mmigr at which galaxies on average came to be pas- 
sive (Fig. [5]). The solid line in Fig. [5] is a by-eye lin- 
ear fit to the data: \og{Mmigr/ Mq) = 10.30 + 0.55^, 
i.e. the population migration progresses from mas- 
sive systems at high z to less massive systems at 
lower z ('population downsizing'). Is there evidence 
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Fig. 5. Typical mass at which galaxies migrate on aver- 
age to the passive sequence at a given z (crosses, see Fig. 
m. The solid line is a eye-fit to the crosses. The dashed 
lines correspond to the mass Mcross at which the late- 
and early- type GSMFs crosses, where different mor- 
phological, photometrical and SSFR criteria were used 
to separate the sample into late- and early-type galaxies 
(Pozzetti et al. 2010). See text for more details. 

from direct observations of a migration mass to the 
passive population as a function of z? As men- 
tioned in §§2.3, observations allow now to deter- 
mine the GSMF at different epochs decomposed 
into early- (red) and late-type (blue) galaxies. For 
example, based on the zCOSMOS survey, Pozzetti et 
al. (2010) determined the mass at which the early- 
and late-type GSMFs cross (different estimators 
for these two populations are used) from z 1 to 
z ss 0.2. Such a crossing mass, Mcross, is interpreted 
namely as the typical mass of late-type galaxies mi- 
grating to early- type ones (see also Bell et al. 2007). 
Local estimates of Mcross by Bell et al. (2003; filled 
triangle) and Baldry et al. (2004; filled square) are 
also plotted, as well as the law inferred by Drory & 
Alvarez (2008) for the mass above which the SFR 
as a function of A/g begins to drop exponentially, 
Mqucnch/Me = 10^°-^^{1 + z)2 i. It is remarkable 
the agreement between direct estimates of different 
definitions of the characteristic mass above which 
the passive (red) population of galaxies dominates 
in number density (sec also Bundy et al. 2006; Hop- 
kins ct al. 2007; Vergani ct al. 2008) and the semi- 
empirical inferences of FAIO. 

Migration rates and the quasar connection.- An- 
other prediction that can be made with the semi- 
empirical approach is the rate in comoving number 
density of galaxies migrating from active to passive 
ones as a function of z. In FAIO it was found that 
this rate, 0g(A/migr), up to z = 1 scatters around 
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Fig. 6. The Mmigr mass and its abundance translated to a 
quasar luminosity, Lq, and abundance at three redshifts 
{z = 0.5, 1, and 2, filled stars). The obtained results 
are showed over the luminosity function plots of quasars 
at the three z bins indicated in the panel (Groom et al. 
2004; Croton 2009; see text for details). The masses and 
abundances of galaxies in the migration process roughly 
coincide with the characteristic Lq of the quasar lumi- 
nosity functions at the same redshifts. 

(1.0 — 5.5) X 10"'' gal/Gyr/Mpc-^ without any clear 
trend with z. These rates arc in agreement with the 
estimated passive (red) population growth rate per 
unit of comoving volume determined by Pozzetti et 
al. (2010) for a redshift bin centered at z = 0.34. 

If the typical time of galaxy transition (or 
quenching) is tq , then the abundance of galaxies per 
unit of comoving Mpc in the process of migration, i.e. 
those of mass Afmigr, is </'g,migr ^ '/'g(A^migr) xig- Let 
assume tg — 100 Myr and let estimate the typical 
luminosity of a bright quasar, Lq{z) associated to 
the stellar mass M^ig^iz) or its corresponding halo 
mass M\i{z) (see e.g., Croton 2009) by making some 
assumptions {Lq = ryiEdd, ^Edd = C x Mbh, Mbh 
is the supermassive black hole mass) and by using 
the observable correlations of A/bh with spheroid ve- 
locity dispersion or stellar mass. The result is very 
encouraging (Fig. [5]): at each z (up to z ~ 3), 
both L(5(Mmigr) and its abundance, 0g,migr, match 
the characteristic luminosity of the two-power law 
QSO luminosity function and its abundance as de- 
termined in Croon et al. (2004) and Roberts et al. 
(2006). This suggests that the transition from active 
to passive regime (quenching) in massive galaxies is 
associated to the QSO active phase. 

3.2. Implications 

In the semi-empirical approach described above, 
the observational input is the total GSM F measured 



at different z, while the information provided by the 
ACDM theory allowed to assign halos to galaxies 
and "individualize" average galaxy Mg tracks as a 
function of mass. By using the inferred Mg tracks 
(GHETs), the SSFRs as a function of and z, the 
evolution of the migration or quenching mass Afmigr, 
and the rate per unit of comoving volume of galax- 
ies transiting from active to passive ones were pre- 
dicted. The fact that these predictions are in reason- 
able agreement with direct observational measures 
at different z implies that both the connection be- 
tween galaxies and ACDM halos, and the predicted 
halo MAHs are consistent with observations. 

The semi-empirical results imply in a natural way 
the bimodal distribution of galaxies into active (blue) 
and passive (red) ones, with a bimodality mass scale 
that shifts downward as time progress (more galaxies 
become on average passive) , this mass scale being at 
z ~ Ms = 2 — 3 X 10^" Mq, in agreement with direct 
observations (e.g.. Bell et al. 2003; KaufFmann et al. 
2003; Blanton & Moustakas 2009). Therefore, mass 
seems to be the main driver of color, SFH, and SF 
quenching, though environment is expected to play 
also a role. Observational studies show indeed that 
color and SFR are moderately modulated by envi- 
ronment for intermediate-mass galaxies, which as- 
sembled more recently, in epochs when larger-scale 
structures collapse, affecting this those galaxies liv- 
ing in such overdense regions (Cucciati et al. 2010; 
Peng et al. 2010). The downsizing of the quenching 
mass has been found also to be slightly accelerated 
in overdense regions (Bundy et al. 2006). 

The inference of SSFR (as well as Mmigr) as a 
function of z was under the hypothesis that the Mg 
growth is driven only by local SF, i.e. the growth 
by accretion of stellar systems (dry mergers) was ex- 
cluded. The rough agreement between predictions 
and direct measurements implies that this hypoth- 
esis is reasonable or, said in other words, that the 
former growth channel dominates. 

4. CHALLENGES AND CONCLUDING 
REMARKS 

The semi-empirical results imply that galaxy for- 
mation inside the growing A CDM halos is consistent 
in general with observations. However, these results 
show that the average stellar mass assembly histo- 
ries of galaxies shift from those of their halos in a 
peculiar way and depending on mass (Fig. [3]). The 
complex astrophysical processes and environmental 
effects involved in galaxy evolution should explain 
this. Among the key challenges, we highlight and 
discuss the following ones: 
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4.1. Why the Ms assembly of low-mass galaxies is 
systematically delayed with respect to the 
assembly of their halos ? 

The SAMs introduce (extremely efficient) SN- 
driven galaxy outflows for reproducing the shallow 
faint-end of the luminosity function or, equivalently, 
the high inclination of Mg/Mh at lower masses (Fig. 
12). Can galaxy outflows explain also the inefficient 
galaxy SF in the past as smaller are the halos? 
By means of self-consistent disk-galaxy evolutionary 
models, Firmani, Avila-Reese & Rodriguez (2010) 
have shown that SN-driven outflows, as well as the 
local SF and ISM feedback processes, may deviate 
the SFR history of disk galaxies from the associated 
halo mass aggregation rate history but not enough 
as to reproduce the semi-empirical and empirical re- 
sults. In Figs. [21 [3l and |4] two evolutionary mod- 
els of low-mass disk galaxies are plotted (dotted-line 
curves). The model tracks differ significantly from 
the GHETs. As seen in Fig. [3J the model Mg tracks 
follow closely those of their halos at late epochs, in 
disagreement with the systematical shift as the mass 
is smaller observed for the corresponding GHETs. 
Although the low Ms/ Mh ratios were obtained in the 
models at z = (by assuming very high SN kinetic- 
energy injection to the outflow, see also e.g., Dutton 
& van den Bosch 2009), at higher z, model predic- 
tion are far from the semi-empirical inferences (Fig. 
[2]). The disagreement between models and observa- 
tions is also evident by comparing the SSFRs: the 
former predict SSFRs for small galaxies much lower 
than the observed ones, specially toward lower z. 

Firmani et al. (2010) explored the possibility of 
later re-accretion of the ejected gas (e.g., Oppen- 
heimcr et al. 2010). For reasonable schemes of gas 
re-accretion as a function of Mh, they found that the 
SSFR of galaxies increases but in the opposite direc- 
tion of the downsizing trend: the increase is large for 
massive galaxies and small for the less massive ones. 
From the side of SAMs, it was found that the popu- 
lation of small galaxies (both central and satellites) 
is too old, red, and passive as compared with obser- 
vations (e.g., Somerville et al. 2008; Fontanot et al. 
2009; Santini et al. 2009; Liu et al. 2010), issues 
that certainly are related to the mentioned-above 
problem of the Ms buildup of ACDM-based sub-M* 
model galaxies. Current N-body + hydrodynam- 
ics cosmological simulations (where feedback-driven 
outflows are allowed) of central low-mass galaxies 
(0.2<Ms/10^Mq<30 at z = 0) show similar prob- 
lems (Colm et al. 2010; Avila-Reese et al. 2011): 
(1) the SSFRs are 5-10 times lower than the average 
of observational determinations at 2 ~ 0, an incon- 



sistency that apparently remains even at z ^ 1 — 1.5, 
though less drastic; (2) the Ms/ My, ratios are ^ 5-10 
times larger than observational inferences at z ~ 
and this difference increases at higher z. 

Unless the current observational inferences of 
A/g, SFR, and their distributions at different z are 
dominated by strong systematic effects and selec- 
tion biases (see for discussions BCWIO; Firmani et 
al. 2010; Avila-Reese et al. 2011), the problems 
highlighted above pose a sharp challenge to current 
ACDM-based models and simulations of low-mass 
galaxies. Its solution requires that galaxies smaller 
than Ms{z = 0)^ 2 X 10^" Mq should have signif- 
icantly delayed their Ms assembly with respect to 
the assembly of their corresponding ACDM halos, 
besides the smaller the galaxy, the longer should be 
such a delay. This is in line with the staged galaxy 
formation picture proposed in iNoeske et al~l pOOTi ). 



4.2. Why the more massive the galaxy, the earlier 
assembled most of its Ms ? 

Though at first glance contradictory to the 
ACDM picture, this manifestation of downsizing has, 
at least partially, its natural roots namely in the hi- 
erarchical clustering as sembly of halos and t heir pro- 
genitors distributions ( Neistein et al.l [20061 see also 
Guo & White 2008; Li, Mo & Gao 2008; Keres et al. 
2009). Besides, red massive galaxies are expected to 
have been formed in early collapsed massive halos- 
associated to high peaks, therefore highly clustered- 
that afterwards become part of groups and clusters 
of galaxies, leaving truncated therefore the early effi- 
cient mass growth of the galaxies associated to these 
halos. The measured correlation function of lumi- 
nous red galaxies is indeed very high (e.g., Li et al. 
2006). On the other hand, massive galaxies typi- 
cally hosted in the past AGNs. The strong feedback 
of the AGN may help to stop gas accretion, trun- 
cating further the galaxy stellar growth and giving 
rise to short er formatio n time-scal es for the mas- 
sive galaxies (Bower et al.i ,2006; .Croton et al. 2006t 
De Lucia et al.ll2006l ). Although all these effects may 
explain the archeological downsizing, several ques- 
tions remain yet open and subject to observational 
testing, among them, whether AGN-feedback is as 
efficient as required for quenching galaxies. 

4.3. Why galaxies transit from active to passive Ms 
growth regimes, the more massive earlier than 
the less massive ones? What does quench the 
SFR of galaxies? 

These questions are related to the previous one 
but make emphasis on the continuous and mono- 
tonic decreasing of Afniigr with time. An important 
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question to take into account is the possible link 
between SFR quenching and morphological trans- 
formation; after all, active (blue) and passive (red) 
galaxies use to be identified with disk- and spheroid- 
dominated systems, respectively. As numerical simu- 
lations show, the main driver of morphological evolu- 
tion seem to be mergers (e.g., Scannapieco & Tissera 
2003; Cox et al. 2006 and more references therein). 
Mergers can also quench SF through the associated 
explosive quasar or starburst phase that heats or 
drive out cold gas in the spheroidal merger rem- 
nant. Thus, the scenario where ACDM-based merg- 
ers drive the formation and quenching of red, early- 
type galaxies passing by a quasar and/or starburst 
phase looks promising (Kauffmann & Haehnelt 2000; 
Hopkins ct al. 2008, and more references therein). 
The connection showed in §§3.1 between Minigr(z) 
and the characteristic luminosity of the quasar lumi- 
nosity function at each z supports this scenario, at 
least for the formation of red galaxies more massive 
than Mmigr (M^igr « 2 10^° M© at z = 0). 

However, recent observations suggest that the sit- 
uation is more complicated, specially for the origin 
of passive (red) galaxies of smaller masses: nearly 
50% of the red-sequence COSMOS galaxies (z<l) 
have disk-like morphologies, being bulge-dominated 
(Bundy et al. 2010), though above ~ 10" M©, this 
fraction is smaller and decreases with time. The au- 
thors suggest that passive disks may be a common 
phase of galaxies migrating to the red sequence, and 
once formed, their transformation into spheroidals 
is moderately fast (1-3 Gyr), driven likely by minor 
mergers. Passive disks could be the result of disk re- 
growth after a gaseous (wet) major merger (Springel 
& Hernquist 2005; Governato et al. 2009), where 
the formed bulge may host an AGN able to quench 
further SF or the bulge can act as a suppressor of 
disk instabilities and SF (Martig et al. 2009). In 
Bundy et al. (2010) are explored several alternatives 
in the light of their observational results; they con- 
clude that likely some combination of processes and 
environmental effects may be operating simultane- 
ously in the formation of the red-sequence galaxies. 

4.4. What are we lacking? 

In case the current observational picture is con- 
firmed with more and deeper observational studies, 
the above-mentioned issues are likely related to our 
lack of understanding of several complex astrophysi- 
cal and environmental processes of galaxy evolution 
in interaction with the dark mould. A key ingredient 
not yet well studied is the physics of the intergalactic 
medium (IGM), taking into account the generation 



and dissipation of turbulence. The trapping and in- 
fall to the center of halos of the IGM, as well as the 
feedback from galaxies (mainly due SNe and AGN) 
with this medium, are relevant phases of galaxy as- 
sembly that crucially depend on the physical state of 
the IGM, especially at early epochs. The SF process 
is another poorly understood ingredient, especially 
in the regime of mergers and strong gas shocks or 
in subcritical conditions (e.g., low gas surface densi- 
ties) . If the above-mentioned issues persist after the 
" bright" side of galaxy formation is well understood 
in the light of observations, then we should have to 
think about possible modifications to the underlying 
cosmological background, the hierarchical clustering 
ACDM scenario. 

V.A-R. thanks the organizers for the invitation. 
We thank PAPIIT-UNAM grant INI 14509 for par- 
tial funding. 
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